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We present a semidefinite program optimization approach to quantum error correction that yields codes and 
recovery procedures that are robust against significant variations in the noise channel. Our approach allows 
us to optimize the encoding, recovery, or both, and is amenable to approximations that significantly improve 
computational cost while retaining fidelity. We illustrate our theory numerically for optimized 5-qubit codes, 
using the standard [5,1,3] code as a benchmark. Our optimized encoding and recovery yields fidelities that are 
uniformly higher by 1-2 orders of magnitude against random unitary weight-2 errors compared to the [5,1,3] 
code with standard recovery. We observe similar improvement for a 4-qubit decoherence-free subspace code 
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Introduction. — Quantum error correction (QEC) is essen- 
tial for the scale-up of quantum information devices. A the- 
ory of quantum error correcting codes has been developed, in 
analogy to classical coding for noisy channels [1-4]. This the- 
ory allows one to find perfect-fidelity encoding and recovery 
procedures for a wide class of noise channels. While it has 
led to many breakthroughs, this approach has two important 
disadvantages: (1) Robustness: Perfect QEC schemes, as well 
those produced by optimization tuned to specific errors, are 
often not robust to even small changes in the noise channel. 
(2) Cost: The encoding and recovery effort in perfect QEC 
typically grows exponentially with the number of errors in the 
noise channel. Here we present an optimization approach to 
QEC that addresses both these problems. Regarding robust- 
ness, we develop an approach that incorporates specific mod- 
els of noise channel uncertainty, resulting in highly robust er- 
ror correction. Thus entire classes of noise channels which 
do not satisfy the standard assumptions for perfect correction 
[1-4], can be tailored with optimized encoding and/or recov- 
ery. Regarding cost, if the resulting robust fidelity levels are 
sufficiently high, then no further increases in codespace di- 
mension and/or levels of concatenation are necessary. This 
assessment, which is critical to any specific implementation, 
is not knowable without performing the robust optimization. 

Relation to prior work. — An optimization approach to 
QEC was reported in a number of recent papers [5-7]. In 
these works error correction design was posed as an optimiza- 
tion problem to directly maximize fidelity, with the design 
variables being the process matrices associated with the en- 
coding and/or recovery channels. Here we present an indi- 
rect approach to fidelity maximization based on minimizing 
the error between the actual channel and the desired channel. 
Both the direct and indirect approaches lead naturally to bi- 
convex optimization problems, specifically, two semidefinite 
programs (SDPs) [8] which can be iterated between recovery 
and encoding. For a given encoding the problem is convex 
in the recovery. For a given recovery, the problem is convex 
in the encoding. However, the indirect approach is in general 
more efficient computationaly then the direct approach and, as 
we will show, has the added advantage of incorporating an ap- 



proximation method which provides a considerable reduction 
in computational cost with only a minimal fidelity loss. 

Organization. — After reviewing the standard error correc- 
tion model and defining performance measures, we state the 
direct and indirect fidelity optimization problems, then show 
how to add robustness measures, and describe methods to 
solve the (robust) indirect problem. We present a range of ex- 
amples illustrating the robustness of our codes against increas- 
ingly more challenging noise channels, and conclude with a 
discussion and summary of computational cost. 

Noise and error correction model. — Subject to standard 
assumptions, the dynamics of any open quantum system can 
be described in terms of a completely-positive (CP) map: 
ps — > ~^2iAipsA\, a result known as the Kraus Operator 
Sum Representation (OSR) [4]. Here ps is the initial sys- 
tem density matrix and the Ai are called operation elements, 
and satisfy ^ A\Ai = Is (identity). The standard error cor- 
rection procedure involves CP encoding (C), error (£), and 
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recovery (1Z) maps (or channels): ps — > pc —> °c - * Ps, 
i.e., using the OSR: p s = Y, r ,e,ci R rEeC c )p s {R r E e C c y 
(see Ref. [9] for a relaxation of the CP map condition). The 
encoding {C c } c fr 1 and recovery {i? r } r= f 1 operation elements 
are rectangular matrices, respectively n c x n s and n s x n c , 
since they map between the system Hilbert space (of dimen- 
sion ns) and the system+ancillae Hilbert space (of dimen- 
sion nc). The error {E e } e J[ operation elements are square 
(n c x n c ) matrices, and represent the effects of decoher- 
ence and noise. As in [3], we will restrict attention to uni- 
tary encoding: C has only a single OSR element, the n c x n s 
encoding matrix C, whose n s columns are the orthonormal 
codewords (C^C = Is) with n c = n s n CA , n CA being 
the dimension of the encoding ancilla space. It follows from 
Ps = ^CA,Ra[Ur(\Q){Q\ra®vc)U r ], where U R is the uni- 
tary recovery acting on the noisy encoded state ac and the 
recovery ancillae |0).ra> that m R = n CA n RA , with n RA the 
dimension of the recovery ancillae space. 

Performance measures. — Our error correction objective is 
to design the encoding C and recovery 1Z so that, for a given 
£ , the map ps — ► ps is as close as possible to a desired 
n s x n s unitary Ls- A common measure of performance 
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is the average fidelity between the channel 1ZSC and the ideal 
L S : /av g = 4E r ,e \Tr LlR r E e C\ 2 . / avg = 1 if and 
only if there are constants a re such that [3, 4]: R r E e C = 
a re Ls, Y2r e \ a re\ 2 = 1- This suggests the indirect mea- 
sure of fidelity, the "distance-like" error (using the Frobenius 
norm, ||X||| = Tr Jf tjQ, 
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the n c x n c m E rectangular error system matrix E 



[Ei ■■■ E mE ), and R the n c n RA 



recovery matrix 



obtained by stacking the n CA n RA matrices R r . Hence, we 
have ||a||f ro = Tr a^a = E r ,eM 2 = 1, and flt R = 
Y, r R\R r = I C - Note that U R = [R\W], where U R is the 
unitary recovery and W is n c n RA x {n c n RA — n c ). 

As / avg and d- in( i are explicitly dependent on the channel 
elements, they are convenient for optimization. Consider then 
the following optimization problems. 



Direct Fidelity Maximization 
maximize / avg = 4^ J^ r e |Tr i^i? r £ , e C| 2 
subject to R j <R = Ic, C^C = I s 
Indirect Fidelity Maximization 
minimize di n d — \\RE(Ie ® C) — a® Ls 
subject to R^R = I C , C^C = I S , 
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(2) 



(3) 



Both are non-convex optimization problems for which local 
solutions can be found from a bi-convex iteration. The direct 
problem was addressed in [5-7]. We now discuss methods to 
obtain local solutions to the indirect problem. 

Robust error correction. — A major limitation of the stan- 
dard procedure of modeling the error channel as fixed, i.e., 
in terms of given operation elements {E e }, is that this does 
not account for uncertainty in knowledge of the channel, and 
in most cases will hence be too conservative. For example: 
different runs of a tomography experiment can yield different 
error channels {£p}p =1 ', an OSRmodel£ (p) could depend on 
an uncertain parameter p; a physical model of the error chan- 
nel might be generated by a system-bath Hamiltonian H (6) 
dependent upon an uncertain set of parameters 9. Whatever 
the source or sources, not accounting for model uncertainties 
typically leads to non-robust error correction, in the sense that 
a small change in the error model can lead to poor perfor- 
mance of the error correction procedure. One way to account 
for uncertainties in terms of an OSR is to take a sample from 
the set, say, {£(pp)}p =1 or {H(9p)}p =v In the latter case, 
tracing out the bath states will result in a set of error system 
matrices {Ep}p_ v To handle this, the objectives in (2)-(3) 
need to be modified. Two possibilities are the worst-case and 
average-case. For the worst-case, these objectives can be re- 
placed by optimizing over all Ep. In the average-case, the 
objectives can be equivalently expressed in the same form but 
with / avg and d m d replaced by their average over all Ep. This 
is equivalent to replacing the error system matrix elements by 

E e I ^ Ep^e/Vi. 



Indirect fidelity maximization. — Using the constraints in 
(3) gives the distance measure (1) as, 
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+ Tr E(I E ® CC ] )E^ - 2Rc Tr RE^ <g> CL\) 



Since only the last term depends on R, minimizing c? inc j over 
R is equivalent to maximizing the last term over R. A singular 
value decomposition of the n c x n c n RA matrix E{oft®CL* s ) 
immediately yields, 



max RcTr RE{cJ ® CL\) = Tr J E(j ® CC^Ei 



R1R=I C 



with the niE x vhe matrix 7 = a^a. The 
optimizing recovery matrix R is given by: 

R=[vi ■ ■ ■ V nG ] [ui ■ ■■ U„ c ] 



(5) 

n C n RA X n C 



(6) 



where the Vi and itj are, respectively, the right and left singular 
vectors of the matrix E(a) ®CL\). Given (C,7), and the fact 
that a need only be chosen so that cvot = 7, the following 
choice for a achieves 7. 
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R is n s n CA x n s n CA (unitary) 



a = V7 
R is n Q m 
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Result (7) implies that R is unitary when the number of en- 
coding ancillas, n CA , is chosen large enough that no recovery 
ancillas are needed, i.e., n RA = 1, and Ur = R. Related re- 
sults about unitarily recoverable codes were obtained in [10]. 
When there are insufficient encoding ancilla, i.e., n CA < m E , 
(8) reveals that additional recovery ancilla are needed so that 
m E = n CA n RA . The result in (7)-(8) does not change if a 
multiplied by a unitary. This unitary freedom is exactly the 
unitary freedom in choosing the OSR operators [4]. Note also 
that the relation n RA n CA = m E may require that R is padded 
with zeroes. 

Optimal recovery. — Given an encoding C, an optimal re- 
covery R can be obtained in two steps. 

Step 1 : solve for 7 which maximizes (5), that is, 



maximize Tr y/E^ <g> CCTjW 
subject to 7 > 0, Tr 7 = 1 



(9) 



Step 2: obtain R from 7 via (6)-(8). 
Since the negative of the objective function in Step 1 satisfies 
the second order condition for convexity, and the constraint is 
a convex set in 7, it follows that (9) is a convex optimization 
problem [8]. 

Approximation to optimal recovery. — If the errors were 
random unitaries, i.e., Ei = ^JplUi, where {pi} are probabil- 
ities and {Ui} are unitaries, then the diagonal elements of the 
matrix 7 would correspond to the probability of the associated 
error [11]. Generalizing to arbitrary channels, we consider the 
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approximation of setting 7 equal to the diagonal matrix with 
diagonal elements, 

lu = \\E t \\lJn c (10) 

each being the average sum-square of the singular values of 
Ei. Since £ is trace-preserving, ^\ la = 1 as required. Us- 
ing this approximation in (7)-(8) to calculate R directly is ob- 
viously very efficient, especially for large dimensions, where 
by comparison solving (9) can be computationally expensive 
even when 7 is constrained to be diagonal. 

Optimal encoding. — Given (R,a), an optimal encoding 
can be found by solving (3) for C. Replacing the non-convex 
equality constraint, C'C — Is, by the convex inequality 
C^C < Is leads to the relaxed convex optimization problem, 

minimize d in d = \\RE(I E ® C) - a ® L s \\% nu 
subject to CtC <I S ' 

By replacing the singular values of the optimal relaxed solu- 
tion to (11) with ones, we obtain a nearby encoding which 
satisfies C'C = Ig. 

Examples. — We now apply the methods developed above 
to the goal of preserving a single qubit (n s = 2) using 4 and 
5-qubit codespaces. We consider noise channels with weight- 
2 and 3 errors, and compare optimal encoding and recovery to 
the performance of the [5, 1, 3] code, which is perfect against 
arbitrary weight- 1 errors [2], and to the 4-qubit DFS code 
(denoted "DFS-4") which is perfect against collective errors 
[12]. We consider an independent errors model, where an er- 
ror on t > qubits has probability P(t) =p t (l — p) n c~ t /Z, 
Z = J27=o where w is the weight. For 5 qubits with 
weight-2 errors there are m E = 2~2t=a (t) = ^ OSR error 
system elements. For weight-3 errors m E = 26. 

Optimized recovery for weight-2, 3 bit-flip errors. — Fig. 1 
shows f avg (R, E, C) vs. bit-flip probability p with weight- 
2 and weight-3 errors. In this figure the code is always the 
standard [5, 1, 3] code and we compare different recoveries. 



The solid lines show the [5,1,3] code performance with stan- 
dard recovery for arbitrary weight-2 or weight-3 errors, which 
as expected, is not good. The curve with circle markers la- 
beled "optimal weight-3 errors" is for an optimal recovery by 
solving (9) for the optimal 7 using the standard [5,1,3] en- 
coding and then finding R from (7)-(8). The improved per- 
formance over the solid lines is entirely due to the optimized 
recovery. Iteration between encoding and recovery [Eqs. (7)- 
(11)] changed neither the optimized recovery nor the [5,1,3] 
encoding. We also see that the optimal performance decreases 
until p = 0.5 and then increases. The average-case recovery 
indicated by the dotted curve with triangle markers falls below 
the optimal forp < 0.5 and then is identical with the optimal 
for p > 0.5. For weight-2 errors, perfect fidelity over the 
entire p-range is achieved by the optimal and average-case re- 
covery (labeled "optimal & avg-case, weight-2 errors") which 
are found by solving either (9) for the optimal 7 or using 
the approximation (10) and then finding R from (7)-(8). The 
perfect performance is easily understood from the fact that 
the [5, 1, 3] code also perfectly corrects weight-2 bit-flip er- 
rors. Our optimization finds the corresponding ideal recovery. 
When this recovery is applied to the weight-3 errors we get 
the dashed curve with star markers labeled "optimal weight-2 
for weight-3 errors". This curve follows the weight-3 optimal 
recovery for p < 0.5 and diverges thereafter, a phenomenon 
similar to what was reported for amplitude-damping errors in 
[5]. In Fig. 1, in all weight-2 cases the recovery matrix R was 
a 2 5 x 2 5 unitary, or easily reducible to that via a singular value 
decomposition, i.e., no recovery ancillae were needed [recall 
Eq. (7)]. For weight-3 errors, m E = 26 requires two recovery 
ancillas, as per Eq. (8) with n RA n CA = m E . 

Weight-2 random unitary errors. — Here, with the excep- 
tion of results marked by "iterated", we always optimized re- 
covery for the DFS-4 or the [5, 1, 3] code. In the "iterated" 
case we start from the DFS-4 or the [5, 1,3] code and op- 
timize the encoding as well. We examined two cases: (1) 
for the DFS-4 code, a single unitary is randomly chosen and 
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FIG. 3: Weight-2 all-random unitary error for 5-qubit codes. 

used for every one of the 15 errors, thus yielding collective 
weight-2 errors; (2) for the [5, 1, 3] code, 15 unitaries are in- 
dependently randomly selected. While the choice of the DFS- 
4 code is somewhat arbitrary, we picked this code for case 
(1) since the errors in this case are similar to the collective 
decoherence model [12]. Figure 2 shows the mean fidelity 
and standard deviation (error bars) from 100 runs using the 
DFS-4 code with recovery being the inverse of the encoding 
(squares), average-case recovery over the whole p-range (tri- 
angles), iterated average-case (diamonds; encoding optimized 
as well), and iterated optimal-case recovery at each p (circles). 
The last two are obtained by iterating (9)-( 11) with the 7- 
approximation (10) to 5 significant digits. The results for the 
case when all 15 unitaries are randomly selected are shown in 
Fig. 3. Here we start from the [5, 1, 3] code and show stan- 
dard recovery (squares) vs optimized results (same legend as 
in Fig. 2). The insets in Figs. 2 and 3 show the fidelity gain 
relative to recovery via decoding in the DFS-4 case, or rela- 
tive to standard recovery ([5, 1, 3] case). In the iterated case, 
where the encoding is also optimized, the gain is uniformly an 
order of magnitude or more. It is important to note that this 
gain is obtained already for small values of p, where standard 
encoding and recovery are designed to perform well. The en- 
coding matrices in the iterated case have support over all basis 
states and and do not appear to yield previously known codes. 

Discussion. — We have presented an optimization approach 
to quantum error correction that yields codes which achieve 
robust performance. One way to interpret the results presented 
in Figs. 1-3 is to consider a scenario where one is faced with 
a noise channel with weight-2 or 3 errors, and can only use 
4 or 5 qubits to encode one. Without optimization perhaps 
the most reasonable choices are the DFS-4 and [5,1,3] codes. 
However, our results show that one can obtain at least an order 
of magnitude higher fidelities by optimizing both the encod- 
ing and recovery procedures, while respecting the constraint 
of 4 or 5 qubits. There are of course standard codes dealing 
with weight-2 and higher errors, but they require significantly 
more qubits. We stress again that knowing that such perfor- 
mance is possible can alleviate the need for unnecessary ad- 



ditional codespace which may be impractical, and can help in 
improving the fault tolerance threshold. 

We note that the optimization approaches presented here 
have differing computational costs. Evaluating this cost de- 
pends on the optimization algorithm and the problem structure 
[8]. The cost can vary greatly if the algorithm is modified for 
the specific problem structure. For general comparison pur- 
poses, all measures of computational complexity clearly will 
depend on the dimensions of the various search spaces. This is 
summarized in Table I which gives the number of optimization 
variables in the SDP optimizations. Clearly, the approximate 
indirect method enjoys superior scaling. 

An intriguing prospect is to integrate the results found here 
within a complete "black-box" error correction scheme, that 
takes quantum state or process tomography as input and it- 
erates until it finds an optimal error correcting encoding and 
recovery. Another important open problem is to extend the 
robust procedures developed here into a fault-tolerant error 
correction scheme. 
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n s n c 



TABLE I: Number of optimization variables (n c = n s n CA ). 
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